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It is icnown that optical-lattice (OL) potentials can stabilize solitons and solitary vortices against 
the critical collapse, generated by the cubic attractive nonlinearity in the 2D geometry. We demon- 
strate that OLs can also stabilize various species of fundamental and vortical solitons against the 
supercritical collapse, driven by the double-attractive cubic-quintic nonlinearity (however, solitons 
remain unstable in the case of the pure quintic nonlinearity). Two types of OLs are considered, pro- 
ducing similar results: the 2D Kronig-Penney "checkerboard", and the sinusoidal potential. Soliton 
families are obtained by means of a variational approximation, and as numerical solutions. The 
stability of the families, which include fundamental and multi-humped solitons, vortices of oblique 
and straight types, vortices built of quadrupoles, and supervortices, strictly obeys the Vakhitov- 
Kolokolov criterion. The model applies to optical media and BEG in "pancake" traps. 

PACS numbers: 03.75.Lm, 05.45.Yv, 42.65.Tg, 42.70.Nq 

I. INTRODUCTION AND THE MODEL 

Formation of multidimensional solitons and solitary vortices (solitons with embedded vorticity) has drawn a great 
deal of attention in studies of nonlinear optics and Bose-Einstein condensation (BEC), see review [1]. (2-1-1)- 
dimensional spatial solitons and quasi-2D spatiotemporal ones were created in crystals with the quadratic (x^^"*) 
nonlinearity [2]. Also reported were spatial solitons, vortices, and dipole-mode states in photorefractive crystals with 
a photo-induced lattice, where the nonlinearity is saturable [3]. However, truly 2D or 3D solitons have not yet been 
observed in media with the generic cubic (x^^'') nonlinearity, the problem being their instability against the collapse 
[31 [S]. Vortex solitons are subject (in the uniform space) to a still stronger azimuthal instability, which occurs even 
in the absence of the collapse [S] . 

A relevant ingredient of both optical and BEC settings which may stabilize multidimensional solitons and vortices is 
an effective periodic potential. In optics, it implies a periodic modulation of the refractive index in the transverse plane, 
as in photonic-crystal fibers), and in BEC it is induced by optical lattices (OLs). In the medium with the self- focusing 
(SF) x^'^-' nonlinearity, the stabilization of fundamental and vortical 2D solitons (those with topological charges S — 
and 1, respectively) under the action of the square-lattice potential was predicted in works [7] and [8] (see also works 
[S]; stable solitons were found too in photonic-crystal-fiber models [TO'). The action of the stabilization mechanism 
can be summarized as follows. In the free 2D space, the x^'^-' nonlinearity supports a family of Townes solitons (TSs), 
which have the same norm, Q — Qts, at all values of the propagation constants, k [5] [TT]. Although there are 
no unstable eigenvalues in the spectrum of small perturbations around Townes solitons, a specific zero eigenvalue 
accounts for the instability against sub-exponentially growing perturbations, which is a manifestation of the critical 
character of the collapse induced by the x^'^'' nonlinearity in the 2D space. As demonstrated by numerical findings 
and the variational approximation (VA), the action of the OL potential with small strength e lifts the degeneracy of 
the Townes-soliton family, stretching the single point, Q = Qts, into interval Qts — < Q < Qts of width AQ ~ e, 
which is filled by stable solitons, that obey the Vakhitov-Kolokolov (VK) criterion, dQ/dk > 0. This criterion is well 
known to be a necessary condition for the stability of solitons in media with SF nonlinearities [51 112j. 

As concerns vortex solitons supported by the square OL, the most compact "crater-shaped" ones, represented by 
a single density peak with a vorticity-induced hole in the center, are unstable [13J. Stable vortices with S = 1 can 
be built as sets of four (or eight) peaks, with the phase shift tt/2 (or 7r/4, respectively) between adjacent ones [TJIH]. 
Stable vortices of higher orders, up to = 6, built of up to 12 peaks, were found too, as well as stable "supervortices", 
i.e., ring chains of 12 (or more) compacts vortices carrying local spins s = 1, with global vorticity S = ±1 imprinted 
onto the entire ring. These patterns were found in models with cubic and saturable SF nonlinearities |15j . Stable 
vortices of the gap-soliton type were also found in the model combining the OL and self-repulsive nonlinearity |16j . A 
general mathematical classification of various localized states with an intrinsic phase structure (in particular, vortices) 
in the 2D model with the square-lattice cos potential and cubic nonlinearity was developed in [17] . 

In the 2D geometry, supercritical collapse is generated by the SF quintic (x*'^^) nonlinearity, which may occur in 
a combination with x^'^'' terms. In optics, the cubic-quintic (CQ) nonlinearity with the SF x'^^ part was predicted 
[H] and recently observed [TOj in aqueous colloids. It was also observed in dye solutions [SU], and recently predicted. 
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through the cascaded mechanism, in two-level media It is relevant to mention that the x*-^-* nonlinearity with 
the self-defocusing sign is observed in various uniform optical media |22j. 

Besides the context of nonlinear optics, the CQ nonlinearity of the SF type appears in the description of BEC with 
attraction between atoms trapped in a "pancake" configuration. Strictly speaking, the reduction of the underlying 3D 
Gross-Pit aevskii equation (GPE) to the 2D form produces a more complex nonpolynomial nonlinearity [23 . However, 
for ID configurations corresponding to a cigar-shaped trap, it was demonstrated that the respective GPE with the 
SF quintic term, which represents a "vestige" of the underlying multi-dimensionality, provides for an appropriate 
description of the soliton dynamics, unless one is interested in the collapse per se |24j . Another straightforward 
interpretation of the quintic term in the GPE is the contribution of three-body collisions, provided that the lossy part 
of this interaction (kicking out atoms from the condensate) may be neglected [25] , 

The OL potential is not necessary for the stability of ID solitons in the SF CQ model, as they are stable in the 
free ID space, despite the possibility of the collapse PB]- On the contrary, in the 2D model with the SF x^'^^ and 
X*^^) nonlinearities, the OL is a crucial factor for the stabilization of solitons against the supercritical collapse. In the 
normalized form, the corresponding model is based on the equation for the mean-field wave function, u{x, y, t) (or the 
local amplitude of an electromagnetic wave, in terms of optics): 



iut + Uxx + Uyy — V{x, y)u + 2|upu -|- 7|u|'*u = 0. (1) 

The derivation of Eq. (jlj from the full GPE yields 7 ~ a1/ , where is the width of the transverse confinement, 
and D the half-period of the OL potential. We will consider two different potentials, viz., the cosinusoidal (cos) 
one, V = — Vo [cos (nx/D) + cos (ny/D)], and the 2D Kronig-Penney (KP) potential, in the form of a "checkerboard" 
composed of cells of size D with potential difference Vq between adjacent ones [T3] (2D states in the checkerboard 
potential combined with the CQ nonlinearity in which the x^^-* term is self-defocusing, i.e., 7 < 0, hence the collapse 
does not occur, were recently investigated in Ref. [13 ). Equation ([ij conserves the norm, 

Q = J J Hx,y)\^dxdy. (2) 

Undoing scalings used in the derivation of Eq. Q from the underlying GPE, one concludes that the actual number 
of atoms in the BEC is ~ (oz/ |as|) Q, where is the scattering length. In terms of nonlinear optics, evolution 
variable t in Eq. ([T]) is the propagation distance, rather than time. 

For the purpose of comparison, we will also consider a modification of Eq. (nj) with pure quintic nonlinearity. 



iut + Uxx + Uyy ~V{x,y)u+ {ul'^u ^ 0. (3) 

Stationary solutions to Eqs. ([T]) and ^ for 2D fundamental (single-humped) solitons were found by means of the 
VA and in a numerical form, as reported below in Section II. Higher-order (multi-humped) nontopological solitons, 
as well as vortices (including those built of quadrupoles, rather than of fundamental solitons) and supervortices, were 
constructed by means of numerical methods, see Sections II and HI, respectively. The stability of all these localized 
patterns was inferred from the VK criterion and verified by direct simulations. To test the stability numerically, 
the initial perturbation was, typically, imposed by the multiplication of a numerically exact stationary solution by 
perturbing factor (1 -I- e), typically with |e| < 0.03 (this was quite sufficient to identify stable and unstable solitons, 
monitoring their evolution in the course of sufficiently long simulations). As a result, a border between stable and 
unstable solutions has been found in each soliton family in the CQ model, while all solitons of Eq. ^ are unstable. 
Approaching the stability border, the width of the stability margin, i.e., the maximum value of perturbation amplitude 
e in the above expression, which does not trigger the onset of instability, shrinks to zero. 



II. FUNDAMENTAL AND MULTI-HUMPED SOLITONS 

Solutions to Eq. (nl) with real chemical potential ~k (in optical models, k is the propagation constant) are sought 
for as u{x, z) = exp [Tkt) U{x, y) , with U {x, y) obeying the stationary equation, 

Uxx + Uyy - V{x, y)U + 2\U\^U + -f\U\*U = kU, (4) 

which is associated with Lagrangian L ^ J J L{x,y)dxdy, whose density is 

L = \Ux\' + \Uyf + [k+V{x,y)] \uf - \uf - (7/3) \Uf . (5) 
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FIG. 1: (Color online) A typical example of numerically found stable fundamental soliton in the model with the KP 
(checkerboard) potential, for £> = 3, Vb = 2, and 7=1. The soliton is pertains to fc = 1.53, its norm being Q = 2.05. 



Fundamental solitons with amplitude A and width W may be approximated by ansatz U {x,y) ~ 
Aexp (— (x^ + y^) / (2VK^)), whose norm (|2|) is Q = nA'^W^, hence the ansatz can be written as 



U{x,y)^ exp , (6) 



W ' \ 2W 

where Q and W may be treated as free variational parameters. The use of the isotropic ansatz is suggested by 
numerically found shapes of the fundamental solitons, which feature an approximate axial symmetry, except for 
"tails" , where the local amplitude of the wave field is very small, see Fig. [l] 

Subsequent derivation of the VA from Lagrangian ^ and ansatz ^ is straightforward for the cos potential. In 
the KP model, the checkerboard potential may be replaced, for this purpose, by its two lowest spatial harmonics 
[13]. Thus, the calculation of the respective effective Lagrangian, L — L{Q,W), leads to the variational equations, 
dL/dQ = dL/dW = 0, which take the form of 



2 



1 + exp 



4L>2 



= 0, 

(7) 



2tt 97r2M^2 g£,2 ^^Pl 4£,2 



The modification of the VA for the quintic-only equation ^ is obvious: in Eqs. ([T]), the terms linear in Q should be 
dropped, and 7=1 should be substituted. 

Families of fundamental solitons, as predicted by the VA, i.e., obtained from a numerical solution of Eqs. (|7]), and 
found from a numerical solution of Eq. Q, that was performed by means of a modification of the relaxation method, 
are represented in Fig. |2]by a set of curves Q{k), for different values of x^^^ coefficient 7. The case of the pure quintic 
nonlinearity, which corresponds to Eq. ([3|, is included too. 

It is seen that solitons in the CQ model do not exist with the norm below a threshold value, Qmin, which is explained 
by the delocalization transition [14 in the region dominated by the cubic nonlinearity (therefore Qmin weakly depends 
on 7). Further, the VK criterion suggests that, in the presence of the cubic term and for any 7 > 0, the solitons 
are stable only up to a point at which dQ/dk changes the signs. In accordance with works [7], the stability-change 
point does not exist in the model with 7 = 0. Note that the norm of stable solitons attains the largest value, which 
is Q = QtS) for 7 = 0; with the growth of 7, the largest norm of the stable solitons decreases, which demonstrates 
the increasing difficulty in the stabilization of the solitons with the transition from the critical (cubic) to supercritical 
(quintic) nonlinearity in the 2D geometry. In accordance with this trend, all solitons in the pure quintic model (the 
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FIG. 2: The norm of the fundamental soHtons versus the propagation constant, k, as predicted by the VA (dotted 
curves) and found from numerical solution of Eq. Q with the cos potential (solid curves) . Fixed parameters are 
Vq — 2 and D = 3. In addition, the dashed-dotted curve displays the same dependence as predicted by the VA for 

the model with the pure quintic nonlinearity, i.e., Eq. ([s]). 

one without the cubic term) are unstable (as per the VK criterion), as the Q{k) curve for the quintic model (the 
dashed-dotted curve in Fig. [2| features solely the negative slope, dQ/dk < 0. Direct simulations exactly corroborate 
all the predictions of the VK criterion. While Fig. |2] displays the results for the model with the cos potential, the 
situation in the KP model is the same, the respective Q{k) curves being very close to those shown in Fig. |2] 

In the CQ model, the decrease of Vq leads to shrinkage of the portions of the Q{k) curves with the positive slope, 
and they disappear at some (Vb)„jj„, leaving only unstable states, with dQ/dk < 0. The respective stability regions 
for the fundamental solitons in the (Vb,Q) plane are displayed in Fig. [s] for the KP model (the situation for its 
counterpart with the cos potential is very similar). The lower stability boundary in this figure is the delocalization 
border, below which solitons do not exist, while the upper border is exactly predicted by the VK criterion (i.e., solitons 
are unstable above it). The evident trend to the shrinkage of the stability region with the decrease of D is explained 
by the exponential smallness of the force of the interaction of a broad soliton with a short-period OL. 

Values of k for all solutions reported in this work belong to the semi-infinite gap in the OL-induced spectrum. 
Solitons can also be found in finite bandgaps; however, as well as in the 2D KP model with the self-repulsive x''^'* 
term [13], all gap solitons turn out to be unstable. 

As shown above, the supercritical collapse imposes an upper bound on the norm of stable fundamental solitons 
[undoing rescalings leading to Eq. ([ij, one can conclude that the number of atoms in the respective matter- wave 
soliton is < 10^]. Stable localized states with a larger norm can be built as multi- humped solitons. Due to the 
symmetry imposed by the OL in two dimensions, the first species of that type following fundamental solitons features 
five peaks, cf. Ref. [13]. Figure |4] displays examples of stable five-peaked solitons found in the KP and cos models, 
for a common value of the chemical potential, k — 1.8 (the solitons are rotated relative to each other by angle 7r/4 
due to the difference in the definition of the periodic potential in the models) . As well as the fundamental solitons, 
families of these solutions feature the stability-change point, separating portions of the respective Q{k) curves with 
dQ/dk > and < (not shown here). Actually, dependences Q{k) and stability regions in the plane of (Vb,Q) for 
the five-peak solitons are similar to those for the fundamental solitons, which are displayed above in Figs. [2] and |3] 
with a difference that Q is larger by a factor of ~ 4. Injection of more norm gives rise to higher-order solitons. For Q 
still larger than the maximum value admitted by the five-peak solitons, their nine-peak counterparts were obtained. 

Besides the solutions built as complexes of in-phase peaks, the model also gives rise to stable dipole, quadrupole, 
and multi-pole localized states. In particular, examples of a quadrupole can be seen below in Fig. [8j as building 
blocks used to compose a new type of vortices {quadrupole vortices). 
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FIG. 3: Stability regions (between the upper and lower borders) for fundamental solitons in the plane of the OL 
strength (Vo) and soliton's norm (Q) in the KP model for different values of OL half-period D. The coefficient in 

front of the quintic term is 7 = 1. 




FIG. 4: (Color online) Typical examples of the five-humped soliton found in models with the KP (b) and 
cosinusoidal (b) potential, respectively. Both solitons pertain to propagation constant A; = 1.8 and differ in the total 
norm: Q = 12.42 for panel (a), and Q = 14.44 for (b). Lattice parameters are the same as those corresponding to 

Fig. [11 



III. VORTEX SOLITONS, QUADRUPOLE VORTICES, AND SUPERVORTICES 

Solitary vortices with topological charge S are found as complex solutions to Eq. Q with the phase circulation of 
2ttS. Compact ("crater-shaped") vortex solitons, with the vorticity nested in a single peak, are unstable (the "crater" 
splits into a set of nonsteady pulses resembling fundamental soltions, with a single one surviving the subsequent 
evolution). Two species of stable vortices with S — 1 have been found in the present models (with the KP and cos 
potentials alike), either one being arranged as a set of four peaks with phase shifts tt/2 between them. Referring to the 
orientation of diagonals connecting the opposite peaks, which compose the vortices, relative to the KP "checkerboard", 
the species may be called oblique and straight, examples of which are displayed together in Fig. |5] While oblique 
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FIG. 5: (Color online) Juxtaposition of contour plots showing density levels, \U {x,y) in straight and oblique 
vortices (solid and dotted curves, severally). The vortices appertain to Vq = 10 and k = 10, with norms 9.79 and 

9.67 for the oblique and straight ones, respectively. 



vortices include a nearly empty site at the center, the straight vortex places its center at a local potential maximum, 
without any vacancy. 

Families of the solitary vortices of both types are presented in Fig. [6]by the respective Q{k) curves, which summarize 
numerically obtained results. The VK criterion is only a necessary condition for the stability of vortices, as it does 
not detect azimuthal instabilities [1]. Nevertheless, systematic simulations have demonstrated that the stability of 
the vortex solitons in the present model precisely obeys the VK criterion, i.e., they are stable up to the point where 
dQ/dk changes its sign, see Fig. |6] This is explained by the fact that the lattice potential is strong enough to suppress 
the azimuthal instability of the vortices, similar to the situation in the cubic model [7 . Note that, as well as in the 
case of fundamental solitons, cf. |2] the stability-change point, dQ/dk — 0, does not exist in the limit of the cubic 
equation, i.e., 7 = 0, and the norm of the solitary vortices attains its maximum just in this limit. 

The model also supports more complex stable vortex structures, such as supervortices, see an example of the oblique 
type in Fig. [7] Although each compact crater-shaped vortex, of which the structure is built, is unstable in isolation, 
the ring formed by four of them, with the global vorticity imprinted onto it, is stable, as verified by direct simulations. 
Q{k) curves for families of the supervortices are similar to those shown in Fig. J6j their stability also precisely obeying 
the VK criterion. In particular, for the same values of parameters as in Fig. (Trthe stability border (which coincides 
with the point where qQ/dt vanishes) is found at fc « 12. At this point, the norm of the supervortex attains its 
maximum, Q 

max 

« 27. 

Families of higher-order vortices with S > I have been found too. A novel type of stable solitary vortices can be 
constructed using, as building blocks, a set of four quadrupoles (rather than simple peaks, cf. Fig. [5]), with the phase 
shift of 7r/2 between adjacent ones, which corresponds to S* = 1. An example of a corresponding quadrupole vortex of 
the straight type is displayed in Fig. [8j These findings will be reported in a detailed form elsewhere. 



IV. CONCLUSION 



We have demonstrated that various species of solitons and vortices can be stabilized by periodic potentials in 
the 2D geometry against the supercritical collapse, driven by the self-attractive cubic-quintic nonlinearity. For the 
Kronig-Penney and cos potentials, soliton families were obtained by means of the variational approximation and 
in a numerical form. The stability of all the families, including fundamental and multi-humped solitons, solitary 
vortices of the oblique and straight types and supervortices precisely obeys the VK criterion. A novel species of 
quadrupole vortices has been demonstrated. The model can be realized in composite (colloidal) optical media, and 
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FIG. 7: (Color online) The distribution of the local density (a) and phase (b) in a stable supervortex supported by 
the KP potential with Vq = 10, D — 2, and 7 = 1, for fc = 9 and Q — 24.26. The global vorticity of the pattern is 
S = +1, while spins of four individual crater-shaped vortices, with centers placed at points {x,y) = (±4, ±4), are 

s = -l. 



in self-attractive BEC in "pancake" -shaped traps. 
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